Limiting fragmentation and meson production at RHIC 
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Abstract 

The PHOBOS's limiting fragmentation etc. three empirical scaling rules for charged multiplicity 
in Au + Au collisions at RBIC are investigated by a hadron and string cascade model LUCIAE. 
Similar studies are performed for the (ft meson exploring its production mechanism via comparing 
with the charged multiplicity. The LUCIAE results for charged multiplicity are compatible with 
PBOBOS observations. Bowever, for the </> meson the three empirical scaling rules are either 
kept only or kept better in the LUCIAE calculations without reduction mechanism of the s 
quark suppression extra introduced for the strangeness in LUCIAE model. These results seem 
indicating a universal production mechanism for charged particle and <ft meson in string fragmenta- 
tion regime. It is discussed that the PBOBOS's empirical scaling rules are model dependent indeed. 
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FIG. 1: Shifted pseudorapidity density distribution of the charged particle per participant pair in 
6% (upper panel) and 35-40% (lower panel) most central Au + Au collisions at - v /s nn =19.6, 130, 
200 GeV. 

The ansatz of limiting fragmentation, i.e. the number of particles produced in high 
energy elementary and/or nucleal collisions by the wounded projectile nucleons should be 
independent of the details of the target, projectile, and beam energy, was introduced very 
early [1]. It was employed later in the fragment production in intermediate energy heavy ion 
collisions. Recently, the BRAHMS and PHOBOS collaborations revealed sequentially the 
evidence of limiting fragmentation in shifted pseudorapidity density distribution of charged 



particle per participant pair in Au + Au collisions at RHIC energies 



Meanwhile, the 



PHOBOS collaboration reported even other two empirical scaling rules: a striking similarity 
between Au+Au and e + e~ collisions in the energy dependence of charged particle production 

n 

and the approximate participant pair scaling of the charged particle production [4] . These 
empirical observations may suggest a universal mechanism of particle production in strong 
interaction system controlled mainly by the available energy, y/s nn , and challenge theorists 
to an explanation. 

On the other hand, there was evidence that the <fi meson production may be distinguish- 



able from the charged particle P, 0]. In addition, the meson is not only a promising 
signature for the QGP formation but also a good probe studying the reaction dynamics at 
early stage [l, 0, O]. It is worthy to investigate above three empirical scaling rules for 
meson exploring its production mechanism via comparison with charged particle. 

To this end, a hadron and string cascade model, LUCIAE [nj], is employed studying 
the limiting fragmentation etc. three empirical scaling rules both for charged particle and 
the (f) meson in Au + Au collisions at RHIC energies in this letter. The three empirical 
scaling rules are reproduced fairly well for charged particle. A discussion is given for the 
meson production mechanism in string fragmentation regime. It is pointed out that since 
the PHOBOS observations rely strongly on the definition and calculation of the number of 
participant nucleons and the PHOBOS's < N part > was extracted based on HIJING jllj], 
the three empirical scaling rules are model dependent indeed. 

The LUCIAE model Q is based on FRITIOF Q- FRITIOF is an incoherent hadron 
multiple scattering and string fragmentation model where the nucleus-nucleus collision is de- 
picted simply as a superposition of nucleon-nucleon collisions. What characterizes LUCIAE 
beyond FRITIOF are the follows: First of all, the rescattering among the participant and 
spectator nucleons and the produced particles from string fragmentation are generally taken 
into account lJ. Secondly, the collective effect in gluon emission of string is considered by 
firecracker model ^| . Thirdly, a phenomenological mechanism for the reduction of s quark 



suppression in the string fragmentation is introduced jl^j ]. 



Fourth, the nuclear shadowing 



effect [la] is taken into account. As the (f> meson does not interact strongly with hadronic 



TABLE I: The charged multiplicity (\u\ <4.7) in 6% most central Au + Au collisions at - v /s nn =19.6, 
130, 200 GeV. 



y/s^. (GeV) 


PHOBOS 1 


LUCIAE (d.) 2 


LUCIAE (m.) 3 


a 


b 


STAR (dN^/dy) 


LUCIAE (dN^/dy) 


19.6 


1670±100 


1466 


1572 


0.5 


0.38 






130 


4020±200 


4779 


4191 


0.05 


1.16 


5.73±1.06 4 


5.08 (6.51 5 ) 


200 


4810±240 


5949 


4964 


0.02 


1.46 


7.20±0.40 6 


8.02 (10.2) 



1. from [3j 

3. modified a and b parameters 
5. with default a and b parameters 



2. default a and b parameters 
4. -0.5< y <0.5 taken from [2PJ 
6. -0.5< y <0.5 taken from [9] 
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FIG. 2: Shifted pseudorapidity density distribution of <f> meson per participant pair in 6% most 
center Au + Au collisions at yjs nn =19.6, 130, 200 GeV from LUCIAE calculations with (upper 
panel) and without (lower panel) reduction of s quark suppression. 
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we neglect it in this work for the moment. 



and 6 meson 



model reproduced fairly well the experimental data of the charged multi 



2l| and the enhanced production of strange baryons (A, H, and O) 



15 



2l| in nucleus-nucleus collisions at SPS energy. However, the LUCIAE model 



overestimates the charged multiplicity, for instance, nearly a factor of 1.2 in Au + Au col- 
lisions at RHIC energy because some energy dependent physics may not well represent in 
LUCIAE model. 

Recently an energy dependent modification of the jet fragmentation function accounting 
for the energy dependence of parton energy loss was proposed in pQCD studies of eA and AA 



collisions \2c 



2^,|2^|. On the other hand, in 



26 



271 128| the string fragmentation function 



in transport model has been considered to be modified in the dense string environment at 



271 . 12811 we assume 



early stage of the relativistic nucleus-nucleus collisions. Following 
that the default a and b parameters in LUND string fragmentation function [29] are suitable 
for Pb + Pb collisions at SPS energies. However, for Au + Au collisions at v /s nn =19.6, 56, 



A+A at SPS and RHIC energies 
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FIG. 3: Charged multiplicity (<j> meson yield) per participant pair in Pb+Pb and Au + Au collisions 
at SPS and RHIC energies as a function of ^/ s nn . 

130, 200 GeV we first adjust roughly the a and b parameters to the experimental charged 
multiplicity |3j. The fitted a and b parameters are then employed to investigate the three 
empirical scaling rules both for the charged particle and meson in above Au + Au collisions. 

In the LUCIAE calculations for Pb + Pb and Au + Au collisions at SPS and RHIC 
energies the centrality, rapidity (pseudorapidity) , and p t cuts are the same as that in the 
corresponding experiments, respectively. As an example, in Tab. 1 the data of charged 
multiplicity (\rj\ <4.7) in 6% most central Au + Au collisions at v /s nn =19.6, 130, and 200 
GeV (taken from jjjj]) are given together with the corresponding LUCIAE results and the 
fitted a and b values. In the LUCIAE results the "d." and "m." in bracket refer to the 
LUCIAE calculation with default (a=0.3 and 6=0,58) and modified a and b parameters, 
respectively. The STAR data of <j) meson rapidity density (-0.5< y <0.5) in Au+Au collisions 
at \/s nn =130 and 200 GeV (taken from (3 and U and the corresponding LUCIAE results 
are given in last two columns. 

In Fig. ^the experimental limiting fragmentation behavior of charged particle in Au + Au 
collisions at Js nn = 19.6, 130, and 200 GeV [3] is compared with the corresponding LUCIAE 
results. The upper and lower panels in Fig. [T] are, respectively, for 6% and 35-40% most 
central collisions. In Fig. [T]the circles, triangles, and squares are, respectively, for ^/s nn =200, 
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FIG. 4: The < N par t > scaling of total charged multiplicity (0 meson yield) per participant pair 
in Au + Au collisions at ^/s rm =200 GeV. 

130, and 19.6 GeV and the full and open labels are for PHOBOS data and LUCIAE results, 
respectively. The shifted pseudorapidity rj is equal to r]-y beam where ybeam refers to the 
beam rapidity (assuming similar value for the pseudorapidity and rapidity variables jsj]). 
One sees in Fig ^ that the LUCIAE model reproduces fairly well the experimental limiting 
fragmentation of charged particle in 6% most central collisions. However, in 35-40% central 
collisions the LUCIAE results are systematically lower than the PHOBOS observations in the 
region rj' > -1.0 . That may attribute to the fact that the discrepancy in < N part > among 
models (PHOBOS's < N p nrt > was extracted based on HIJING model [ll | ) is increased with 
the decrease of centrality j3l|. The PHOBOS's < N part > is visibly lower than LUCIAE in 
35-40% central collisions indeed. 

Given in upper panel of Fig. |21 is the shifted pseudorapidity density distributions per 
participant pair of </> meson in 6% most central Au+Au collisions at, respectively, A /s^=200, 
130, and 19.6 GeV (circles, triangles, and squares) from LUCIAE. Since in LUCIAE model 
a mechanism for the reduction of s quark suppression in the string fragmentation was extra 
introduced for strangeness thus the corresponding LUCIAE results without this mechanism 
is plotted in lower panel. Comparing the lower panel with upper panel we see that the 
limiting fragmentation is kept in a wider r\ region in the LUCIAE calculations without 



reduction of s quark suppression than with ones. 

The data of total charged multiplicity per participant pair in 7%, 7%, and 5% most 
central Pb + Pb collisions at 40, 80, and 158 A GeV (open circles, taken from |4j) and in 
6% most central Au + Au collisions at v /s nn =19.6, 56, 130, and 200 GeV (full circles, taken 
from |4|) are plotted as function of ^/s nn in Fig. El In Fig. El the full and open triangles 
are the corresponding results from LUCIAE calculation with and without the reduction of 
s quark suppression, respectively. The LUCIAE results for meson are given by full and 
open squares (scaled by 100) for, respectively, with and without the reduction of s quark 
suppression. Fig. El turns out that the LUCIAE model fairly reproduces the experimental 
y/s nn dependence of the charged multiplicity per participant pair in Pb + Pb and Au + Au 
collisions at SPS and RHIC energies. For <p meson, only the LUCIAE calculations without 
the reduction of s quark suppression depend on ^s nn in nearly the same way as the charged 
particle. 

In the left panel of Fig. 0] the PHOBOS observation of approximate < N part > scaling 
for total charged multiplicity in Au + Au collisions at v /s nn =200 GeV |4J is given by full 
circle with error bar. The corresponding LUCIAE results with and without reduction of 
s quark suppression are given by full and open triangles. Full and open squares are the 
results of meson from LUCIAE calculations with and without the reduction of s quark 
suppression, respectively. One sees in left panel that the <ft meson yields per participant 
pair from LUCIAE calculations without the reduction of s quark suppression parallel to 
the corresponding charged multiplicity better than ones from LUCIAE calculations with 
the reduction of s quark suppression. The left panel shows also that for charged particle 
although the LUCIAE results are compatible with PHOBOS data within error bar, but the 
LUCIAE results violate < N part > scaling stronger than PHOBOS data. That may attribute 
to the discrepancy in <JV P „, > deration and ealcnlation between PHOBOS and LUCIAE 
model as mentioned in |31|]. To this end, in right panel we compare the PHOBOS data to the 
results (open triangles for < N c h > and open squares for < >) from a calculation where 
< N ch > (< Ntf, >) is from LUCIAE but the < N part > from PHOBOS Q (corresponding to 
the same percentile of the total cross section as in LUCIAE calculation). The full triangles 
and full squares in right panel are the same as that in left panel. We see in right panel 
that using a single value of < Nch > (< N^ >) from LUCIAE but dividing it by < N part > 
from deferent calculation, one from PHOBOS and the other from LUCIAE, the resulted 

n 



< N ch > /(0.5 < N part >) (< > /(0.5 < N part >)) depends on < N part > in different 
way. Thus the above PHOBOS observation is model dependent indeed. 

In summary, we have modified the LUND string fragmentation function that its a and 
b parameters are first adjusted roughly to the experimental data of charged multiplicity in 
Au + Au collisions at RHIC energies jj, Q| . The fitted a and b parameters are then employed 
to study the three empirical scaling rules for both the charged multiplicity and the meson. 
It is turned out that the three empirical observations for the charged multiplicity jjj [^J 
could be fairly reproduced by the LUCIAE model. However, for meson it is either kept 
only or kept better in the LUCIAE calculation without the reduction mechanism of s quark 
suppression. It seems an evidence that the meson production may not distinguish from 
the charged particle in string fragmentation regime, except the mechanism of the reduction 
of s quark suppression extra introduced for the strangeness in LUCIAE model. Because the 
PHOBOS observations rely strongly on the number of participant nucleon and PHOBOS's 

< N part > is extracted based on HIJING 11], the PHOBOS's three empirical scaling rules 
are model dependent. 
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